Magnetic field sensors using large cat states beyond the standard quantum limit 
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Measurement devices could benefit from entangled correlations to yield a measurement sensitiv- 
ity approaching the physical Heisenberg limit. Building upon previous magnetometric work using 
pseudo-entangled spin states in solution-state NMR, we present two conceptual advancements to 
better prepare and interpret the pseudo-entanglement resource as well as the use of a 13-spin cat 
state to measure the local magnetic field with a sensitivity beyond the standard quantum limit. 



Many technologies are looking to quantum mechan- 
ics as a way to dramatically improve upon current ca- 
pabilities [TJ [2]. As examples, interferometry [3j HI [5], 
metrology [6], lithography [7], and information processing 
[5J O GEO] are pursuing quantum information techniques 
to make use of the benefits of highly correlated entangled 
states. Even environments which normally use the pseu- 
dopure approach such as Nuclear Magnetic Resonance 
(NMR) p2 112 EGA EH] can exploit a pseudo-entanglement 
resource for highly sensitive magnetic field measurements 

cm- 

In a locally homogeneous magnetic field, an isolated 
nuclear spin will precess according to its Larmor fre- 
quency which depends only on the nuclear species and 
the magnetic field [16 . Consider the (unnormalised) 
state |0) + |1) in a rotating frame, where the states |0) 
and |1) refer to the parallel and anti-parallel spin eigen- 
states. After a time t such a state would evolve into the 
state |0) + exp(i7f?o^) where the gyromagnetic ratio 
7 is known, so the acquired phase can be used to deduce 
the local magnetic field Bq. A set of N isolated spins 
can therefore serve as microscopic magnetic field sensors 
|16j . with a measurement sensitivity proportional to yN. 
This degree of precision is known as the 'standard quan- 
tum limit' [T7]. 

It is possible to exceed this limit by making use of 
quantum entanglement. If we assemble the N spins in 
the state |0) + |1) = |00...0) + |ll...l), this will evolve, 
after a time t, into the state |0) + exp(iNjB t) |1). This 
evolution allows us to determine Bq with an increased 
sensitivity compared to measuring each spin's evolution 
independently. The degree of sensitivity approaches the 
fundamental Heisenberg uncertainty relation [18], in that 
we are theoretically able to have a sensitivity propor- 
tional to N using N particles. 

We recently reported proof of principle experiments 
exploiting pseudo-entanglement in nuclear spin ensem- 
bles [TS]. In this Letter we grow the size of the cat 
state from 10 to 13 spins and address some limitations 
of the previous approach. Specifically, we incorporate a 
polarisation-priming sequence that more intelligently ex- 
ploits the pseudo-entanglement resource and additionally 
disconnect the centre spin during field measurement to 
simplify the field estimation. 



To quickly generate large states such as |0) + |1), re- 
ferred to as 'NOON' g] or 'cat' [M 1201 123 H2] states, 
we chose natural abundance solution-state tetramethyl- 
silane (TMS) as our sensor compound. We sought a cen- 
tral spin-active nucleus distinct from, and surrounded by, 
many chemically equivalent outer spin-active nuclei, al- 
lowing us to address all peripheral spins globally; this 
highly symmetric configuration means that the pulse se- 
quence complexity is independent of the number of spins, 
in contrast with previous work [T9 l [20" j [21] . Roughly 4.7% 
of the molecules consist of one 29 Si spin surrounded by 
twelve 1 H spins (isotopic labeling could in principle be 
used to increase this proportion, but was not used here). 
Such a molecule is capable of hosting a 13-spin 'cat' state. 

Reading out the thermal state on the centre 29 Si spin 
produces thirteen peaks corresponding to the distribution 
of up and down spins in the nearby hydrogen nuclei. It is 
convenient to assign a number, I, to each of these peaks 
corresponding to their 'lopsidedness', that is i = u — d 
where u and d are the number of intramolecular up and 
down proton spins, respectively. 

The basic sequence consists of a Hadamard gate on the 
centre spin to generate the state |0) + |1), a controlled- 
not (CNOT) gate conditional upon the state of the centre 
spin, followed by an evolution delay t before reversing the 
previous CNOT so that the final phase is mapped onto 
the centre spin for readout. This sequence can either be 
applied to a pseudopure state (corresponding to a single 
peak) or the entire thermal state (all thirteen peaks with 
a Boltzmann distribution of intensities) . The two meth- 
ods are computationally equivalent, but starting from the 
thermal state gives a stronger readout signal. The inter- 
nal peaks will usefully pick up phase proportional to their 
lopsidedness [TS] . The application of the sensor sequence 
to the internal lines generates states that are not quite 
'NOON' states — the term 'Many-Some, Some-Many' (or 
'MSSM') was introduced [TS] to describe these states. 

In this scheme the phase acquired by each peak, iden- 
tified by its lopsidedness £, is given by 

<f>/t = B (^ 7 h + 7si) (1) 

where 7 H = 42.577 MHz/T and 7 Si = -8.465 MHz/T 
are the gyromagnetic ratios of 1 H and 29 Si, respectively. 
As NMR experiments are most conveniently described 
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FIG. 1: Priming the 13-spin MSSM states (created 
using the 12 hydrogen spins (white) and the centre 
29 Si spin (blue) in TMS) is critical for large sensors. 
The three 29 Si spectra show how the polarisation of 
each thermally equilibrated MSSM state (red) can be 
amplified according to its lopsidedness with polarisa- 
tion priming (blue). Each primed MSSM state can 
then used for field measurement, where the very lop- 
sided outermost peaks evolve more rapidly (and are 
hence more sensitive) than the inner peaks. 



in a rotating frame [TB] the observed phase depends not 
directly on the magnetic field but rather on its offset -Bo 
from some assumed nominal value. With this in mind we 
can calculate the phase sensitivity increase of TMS to be 
61.4 over an individual silicon spin and 12.2 over a single 
hydrogen spin. 

The Boltzmann distribution of populations of the 1 H 
spins leads to weak intensities for the outermost MSSM 
lines which are the most sensitive to magnetic field. This 
could be addressed by physically or computationally 
24J manipulating the sensor molecule. We now describe 
a simple approach which uses the quantum resource more 
efficiently than simply averaging many measurements. 

Many techniques for polarisation transfer have been 
developed in NMR systems [25]; these work not by in- 
creasing the polarisation, but instead by transferring po- 
larisation from one part of the density matrix to another 
where it can be more effectively used. One simple exam- 
ple is a CNOT gate applied to a high 7 nucleus controlled 
by a low 7 nucleus, which transfers the population differ- 
ence across a transition of the high 7 nucleus to a low 7 
transition, effectively multiplying the polarisation of the 
insensitive nuclei by 7^, the ratio of the gyromagnetic ra- 
tios of the two nuclei. (In conventional NMR experiments 
this is known as Insensitive Nuclei Enhanced by Polarisa- 
tion Transfer (INEPT) [25, 21)]). In our highly symmetric 
molecular sensor, the amplitude benefits of polarisation 



transfer are even greater than this ratio (7^ = —5.03) 
as the CNOT has the effect of increasing the signal in 
a negative direction for each coupled 'up' 1 H spin (neg- 
ative because 7^ is negative) and positive direction for 
each connected 'down' 1 H spin. Explicitly, a peak of lop- 
sidedness I in a polarisation-primed sensor undergoes an 
amplitude magnification A according to 



A(£) = (1 + 7R £) 



(2) 



This means the outermost lines — the most sensitive 
components of the sensor with the poorest thermal 
populations — are those most amplified by the polarisa- 
tion transfer. The integrated intensities of the outermost 
lines display an approximately 60-fold increase over the 
thermally-polarised measurement as expected. 

Polarisation-priming provides another advantage: by 
swapping the polarisation of sensitive and insensitive nu- 
clei, one only needs to wait for the sensitive nuclei to 
rethermalise. Although the details of relaxation pro- 
cesses can be complicated, in fast-tumbling spin-1/2 sys- 
tems dominated by dipolar relaxation the T\ time is in- 
versely proportional to the square of the gyromagnetic 
ratio [15] , so as the polarisation-priming sequence en- 
hances the polarisation it also decreases the rethermali- 
sation time. A field estimation generated from the basic 
sequence with polarisation priming is shown in Figure [T] 

Instrumental^, there will always be some error asso- 
ciated with imperfect frequency detunings. Equation [T] 
assumes that the frequencies of the 1 H and 29 Si channels 
have been chosen correctly, so that they are precisely on 
resonance with their respective nuclei at the nominal field 
strength. In general: 



<f>/t = B (£ 7H + 7si) + (£Sa + fa) 



(3) 



where <5h and Sgi are the frequency offsets of the nuclei at 
the nominal field. One can mitigate the errors generated 
by these terms by systematically removing them. One of 
these offsets (here assumed to be <5h) can be eliminated 
by shifting the nominal field, but it is only possible to re- 
move both terms if the frequencies are set correctly. An 
imprecisely known rotating frame offset leads to inaccu- 
rate field estimations. This requirement can be removed 
by 'disconnecting' (disentangling) the centre spin during 
the phase acquisition delay. Two methods for achieving 
this are introduced in Figure [2] 

Consider how sequence A acts upon the leftmost line. 
In the pseudopure approximation |13j the leftmost line 
in its thermal state is represented as |0) gi |0)® , where 
® indicates a tensor product. A Hadamard gate followed 
by a CNOT gate conditional upon the silicon nucleus 
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It is here that we can disentangle the central 29 Si spin 
from our large cat state by applying a NOT gate to 29 Si 
in the second term, giving 



io> Si io>r+io)siii>r =iQ) Si (io)i 
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FIG. 2: The two field sensor sequences with all im- 
provements (coloured green) presented in this paper. 
The polarisation-priming component begins the se- 
quence with appropriate delays A = 1/4 J for J the 
spin-spin coupling constant. In Sequence A, the cen- 
tre two CNOT gates are 'disentangling' gates de- 
signed to allow only one nuclear species ( 1 H in TMS) 
to evolve during measurement. In Sequence B, two 
7r pulses separated by half the delay t emulates a dis- 
entangling effect by reversing the phase acquired by 
that nucleus for half of the phase acquisition time. 



so that only the 1 H spins will acquire field-dependent 
phases. 

It might seem that this approach would require 
a multiply-controlled NOT gate (a generalised Toffoli 
gate), but with a pseudopure state this is not required. It 
is only necessary to apply a NOT gate to the second term 
and not to the first term. This can be achieved with a 
modified CNOT gate [13], with the evolution time chosen 
to match the separation between the outermost lines in 
the 29 Si multiplct rather than the conventional coupling 
size. 

This simple approach must be modified to work simul- 
taneously with a general set of MSSM lines. With an odd 
number of X H spins, this can be achieved with a conven- 
tional CNOT gate (see Figure |2j\) , which disentangles 
every MSSM line. This approach does not work for sys- 
tems with an even number of 1 H spins. An alternative, 
simpler method (Figure [2j3) uses echoes to refocus the 
inner 29 Si spin rather than disentangling it, and this can 
be applied to both even and odd systems. Under both 
sequences we read out the acquired phase by applying 
the sequence in reverse, and observing the central 29 Si 
spin. 

To test this approach, we applied a small offset to the 
29 Si channel and implemented a full field estimation with 
the original pulse sequence and with the modified se- 
quence B. As shown in Fig. [3j the field estimation now 
gives slightly different results for different lines in the 
multiplet if the original sequence is used, but these im- 
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FIG. 3: Each peak in the 29 Si NMR spectrum(labeled 
according to lopsidedness) can estimate the local 
magnetic field. The precision of each peak's esti- 
mation scales according to its absolute lopsidedness, 
leading to the most precise estimations at the outer- 
most peaks. The colour amplitude chosen reflects the 
gain in signal from the polarisation-priming sequence; 
a value of 1 on each peak estimation represents its 
thermal binomial amplitude (See Equation[2]). A: The 
field estimation is visibly sensitive to a small nonzero 
detuning (3.5 Hz) on the 29 Si spin if it is not disentan- 
gled during measurement. B: The same 29 Si detuning 
does not distort the field estimation when applying 
disentangling sequence B (see Figure [2^3) 



perfections are removed by the modified sequence. We 
then repeated the phase estimation with a wide range of 
silicon channel dctunings (data not shown) and obtained 
indistinguishable field estimations. Sequences A and B 
were both successfully implemented on our original odd 
spin system, trimethylphosphite (data not shown). 

The most obvious drawback of these new sequences is 
a mild sensitivity decrease from (12jh + lSi)/lH to 12 
times that of a single 1 H spin. Because we remove the 
need to accurately measure frequency offsets, however, 
simplicity (and potentially accuracy) is enhanced. It is 
important to chose the correct disentangling sequence so 
that all MSSM peaks are measured simultaneously, for 
reasons that will now be discussed. 
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We can extract many times more information with a 
single scan by considering all the peaks in the spectrum. 
A single NOON state sensor can easily encounter aliasing 
problems; in effect one already needs to know the approx- 
imate field offset to be certain of the results. A full arse- 
nal of MSSM states provides a mechanism to avoid such 
problems. On an quantum sensor with N outer spins, a 
<j> phase rotation on the outermost peak could only be 
aliased with a rotation of <fi + 2kNir for some integer 
k. In quantum interferometric terms [27 , such a sensor 
simultaneously displays both local and global phase dis- 
tinguishability. Such anti-aliasing effects are a desirable 
property of these highly symmetric sensor molecules. 

For even more sensitive measurement, larger sensors 
can be employed, potentially with iterative (or other) 
geometries. To extract information from the outermost 
peaks of such very large sensors, polarisation amplifica- 
tion methods such as Dynamic Nuclear Polarisation [21] 
can be applied in addition to the methods outlined above. 
All sensors can benefit from both the simplified field esti- 
mation afforded by disentanglement methods, and from 
polarisation-priming the pseudo-entanglement resource 
as introduced in this letter. 

In conclusion, we have generated 13-spin pseudo-'cat' 
states for entanglement-enhanced magnetometry. We 
have proposed and applied innovations to improve the 
stability and resolution of the entanglement resource. 
Errors arising from imperfect knowledge of system vari- 
ables are removed by two different disentanglement meth- 
ods, and the overall weighted enhancement afforded by 
polarisation-priming shows approximately a 60-fold in- 
crease in the most sensitive components of the sensor. 
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(EP/D048559/f). The authors thank Vasileia Filidou 
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Oxford. J.J.L.M. thanks St. John's College, Oxford and 
the Royal Society. 



METHODS 

The sample was a 1:1 by volume solution of tetra- 
mcthylsilane and acetone-d6, degassed using freeze- 
pump-thaw cycles, and flame sealed in a 5 mm Wilmad 
LabGlass NMR tube. All NMR experiments were per- 
formed at a temperature of 20° C on a Varian 1NOVA 
600 spectrometer using a broadband tunable X{H} probe 
with a 2 H lock with a 4-step phase cycle to cancel receiver 
errors, n/2 pulse lengths were approximately 27 fjs on 
the hydrogen channel and 17 fis on the silicon channel. 
The spin-spin coupling 3 JsiH was 6.63 Hz. Measured 29 Si 
relaxation times were T2 = 1.2 s and Ti = 25.4 s, while 
1 H relaxation times were T 2 = 1.6 s and Ti = 8.9 s. 

Quantum logic gates were implemented using standard 
NMR techniques [TJ]. Hadamard gates were applied as 



7r/2_ a 7r 2 . C-NOT gates, equivalent to a controlled-phase 
gate surrounded by Hadamard gates on one channel |14j , 
were implemented as two 1 H tt/2 pulses separated by 
a spin echo of length 1/2 J, where J is the spin-spin 
coupling constant. All Z gates were realised as phase 
shifts in the pulses that followed [TJ]. To reduce off- 
resonance and RF inhomogeneity errors, spin-echoes were 
constructed with two simultaneous tt x pulses at times 
1/8 J and 3/8 J, and all pulses were implemented as si- 
multaneous, equal-duration BB1 composite pulses |29j . 
Implementing such pulses used suitable amplitude ad- 
justments and '0-degree' identity gate pulses where re- 
quired. 

Data was apodised with a Hamming filter and Fourier 
transformed using matNMR [30] version 3.9.59. The 
spectra with no phase accumulation delay was phased, 
and that phase correction was applied to all other spectra 
for consistency. Spectra were then exported to Matlab 
for final processing. 
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